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Locally Cold Flows from Large-Scale Structure 
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ABSTRACT 

We show that the "cold" Hubble flow observed for galaxies around the Milky Way 
does not represent a problem in cosmology but is due to the particular geometry 
and dynamics of our local wall. The behavior of the perturbed Hubble flow around 
the Milky Way is the result of two main factors: at small scales (R < 1 Mpc) the 
inflow is dominated by the gravitational influence of the Milky Way. At large scales 
(R > 1 Mpc) the outflow reflects the expansion of our local wall which "cools down" 
the peculiar velocities. This is an intrinsic property of walls and is independent of 
cosmology. We find the dispersion of the local Hubble flow (1 < R < 3 Mpc) around 
simulated "Milky Way" haloes located at the centre of low-density cosmological walls 
to be <jh ^30 km s _1 , in excellent agreement with observations. The expansion of 
our local wall is also reflected in the value of the measured local Hubble constant. For 
"Milky Way" haloes inside walls, we find super-Hubble flows with h looai ~ 0.77 — 1.13. 
The radius of equilibrium (Ro) depends not only on the mass of the central halo and 
the Hubble expansion but also on the dynamics given by the local LSS geometry. The 
super-Hubble flow inside our local wall has the effect of reducing the radius at which 
the local expansion balances the gravitational influence of the Milky Way. By ignoring 
the dynamical effect of the local wall, the mass of the Milky Way estimated from Ro 
can be underestimated by as much as ~ 30%. 

Key words: Cosmology: large-scale structure of Universe; galaxies: kinematics and 
dynamics, Local Group; methods: data analysis, N-body simulations 



1 INTRODUCTION 



' The mystery of the "Coldness of the Hubble Flow" (CHF) 
has remained largely unsolved despite many attempts to 
explain or reproduce it. Early measurements of the radial 
velocities of nearby galaxies found the per turbed Hubble 
flow dispersion to b e ~. .60 k m s _1 ([Sandage et al.l 

1 19721 : ISandagei Il98fj : iPeeblesi Il988( ). Recen t studies give 
an e v en lower value of an ~ 30 km s -1 (|Ekholm et aU 

l200ll : iKarachentsev et all |2002| , 12001 120091 ) This disper- 
sion is significantly "colder" tha n theoretical expectations 
|Peebleslll98Cl ;ls coccimarr (3 120041 ) and observed pairwise ve- 
locity dispersions which are of the order of 300 — 500 k m s _1 
jZehavi et al.ll2002l : [Landvll2002l : [Hawkins et al.ll2003l h 



1.1 Searching for the origin of the CHF 

The scales involving the CHF suggest that it has a cos- 
mological origin. Two main causes have been proposed 
to explain the CHF: background cosmology (in the form 



of a A term) and local environment. The influence of 
background cosmology was fi rst st udied in N-body sim- 
ulations by iGovernato et all l|l997t ) who failed to detect 
low velocity dispersions in two models: VL = 1 (SCDM) 
and a = 0.3 (OCDM). A later paper (Bar vshev et al.l 
l200ll ) proposed the role of the cosmological constant 
in the CHF under the assumption that the accelerat- 
ing effect of A or a dark energy term results in a de- 
cay in the velocity fluctuations, thereby giving a colder 
local environment |Barvshev et all 200 ll; Teerikorpi et al.1 



20081: iPeirani fc de Freitas P.I 120081 : iMaccio et all booa 



Chernin et al.l l2010h On the other hand, early constrained 



N-body simulations of the local group suggested that the ob- 
served low velocity disp ersion somehow reflects our partic- 
ular local environment llvan de Weygaert fc Hoffman] 120001 ; 
Klvp in et al.1 120 03 V This possibility was also explo r ed by 
Hoff man et al.1 ([20081 ') and iMartinez-Vaquero et aU ([20091 ) 



who analyzed constrained simulations of the local group in 
several cosmologies and found the effect of the cosmologi- 
cal constant on local dynamics to be marginal. While these 
studies reproduced to some degree cold Hubble flows they 
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were not able to identify its origin apart from being an in- 
trinsic property of "local-group" environments. 



2 MIMICKING A FROM LOCAL DYNAMICS 

From the above discussion, a picture begins to appear: in 
order to cool the peculiar velocities, an accelerating A-like 
component is needed, yet there are strong indications that 
this does not come from the background cosmology. This ac- 
celerating factor is somehow related to the particular local 
environment of the Milky Way. One may ask under what 
circumstances one could have a local accelerated expansion 
without a A component. The most straightforward exam- 
ple is that of cosmological voids. Their low density makes 
them underdense island universes and their expansion adds 
to the drag co ming from th e global Hubble expansion as 
Hvoid = 1-2-ffo (lSchaadl200"7h . This effect is independent of 
the background cosmology. It is purely the result of the par- 
ticular dynamics of underdense regions and, as we will see 
in the following section, it is not unique to voids. 



2.1 Local LSS geometry and dynamics 

In order to describe the role of the local geometry in the 
LSS dynamics it is illustrative to take a somewhat differ- 
ent a pproach to the Z el'dovich collapse of an overdense re- 
gion |Zeldovichl Il970l ) and focus on the draining an d ex- 
pansio n of underdense valleys a l ong s imilar lines as 
JHH); llcke fe van de Weygaertl l| 19871 ). We start by gen- 
eralizing the notion of a void as an n-dimensional valley 
with n = 3, 2, 1 corresponding to voids, walls and filaments 
respectively. In this inverse dynamical description, mat- 
ter is drained out of three-dimensional valleys (voids) and 
accumulates in two-dimensional (void-intersecting) walls. 
These two-dimensional valleys expand and lose their matter 
into their surrounding one-dimensional (wall-intersecting) 
filaments. Finally matter flows out of filaments into zero- 
dimensional (filament-intersecting) clusters. 

The dynamics inside a wall can be roughly described 
as a two-dimensional expanding valley. An observer situ- 
ated at the centre of a wall will measure recession veloci- 
ties Vrec of surrounding galaxies as the sum of two compo- 
nents: v rec = + Up, where corresponds to the un- 
perturbed Hubble flow and the line-of-sight peculiar veloc- 



ity Vp — «|J L + i>lss is decomposed here into a non-linear 
component ujJ L coming from small-scale gravitational inter- 
actions and a geometrical component v^ ss originating from 
the particular dynamics of the local LSS geometry where 
the observer is located. In the case of walls, the local expan- 
sion and draining of matter results in v\\ ss being on average 
positive. This extra component in the velocity is reflected in 
super- Hubble flows which further cool down the local veloc- 
ity fluctuations, producing a cold dynamical environment. 



3 ARE WE LIVING IN A WALL? 

The Milky Way is the dominant member of the Local 
Grou p, a system of galaxies forming a relatively flat pan- 
cake ljHartwickll20(5oh . Extending perpendicular to the plane 
of the Milky Way, we find a planar association of galaxies 



up to ~15 Mpc and at its edge the Virgo Cluster and its fil- 
amentary network (|Tullv fc Fisherl 1 1987t h Perpendicular to 
the lo cal wall, we find the local void extending up to ~ 40 
Mpc (|Tullv et al.ll200Sl) . Our surrounding LSS geometry is 
consistent with the Milky Way being inside a wall of radius 
r ~ 10 Mpc. It is interesting to note that the disk of the 
MW and its neighbours within 10 Mpc lie perpendicular to 
the plane of the local wall in exce llent agreement with pre- 
dictions from tidal t orque theory rtHovTell 19491 : IWhitelll98"3 : 

23). 



iNavarro et al1l2004l ; lAragon-Calvo et al J 120071 ) 



4 N-BODY SIMULATION AND DM HALOES 

The analysis presented here is based on a dark matter 
N-body simulation with cosmological parameters Q. m — 
0.3, fi A = 0.7, h = 0.73,a 8 = 0.8 and 200 Mpc/h per side. 
Initial conditions were generated for 512 3 , 256 3 , 128 3 and 
64 3 particles. The simulation was run for all particle res- 
olutions from z = 49 until th e present time , z = using 
the GADGET-2 N-body code (|Springelll2005l ). From the fi- 
nal particle distribution of all resolutions we compute the 
density field inside a cub ic grid of 512 voxels per dimensio n 
using the DTFE method (|Schaap fc van de Weygaertll2000h . 

We identify groups in the 512 3 simulation at z — using 
the Friends of Friends group finder with a linking parameter 
b = 0.2. We impose a particle count limit of 100 correspond- 
ing to a halo mass of M > 5 x 10 M©/i _1 . For each group, 
we compute its mass, radius, peculiar velocity, Vmax , shape, 
etc. We then identify "Milky Way" (MW) haloes in the mass 
range 1 x 10 12 M Q < M halo < 3 x 10 12 M and maximum 
circular velocity t> m „ > 20 km s . 



5 FINDING WALLS IN THE COSMIC WEB 

In order to isolate individual walls from the computed 
density field we performed the hierarchical decomposition 
into voids, walls and filamen ts using the Spine Web method 
l|Arag6n-Calvo et al.l 12010a) extended to the hier archical 
formalism described in lAragon-Calvo et all (|2010bl ). In the 
hierarchical SpineWeb framework, each hierarchical level is 
defined by a physical smoothing scale in the linear regime 
and contains the full network of voids and its complemen- 
tary web of walls and filaments. The hierarchical SpineWeb 
provides a complete parameter-free characterization of the 
geometry of the Cosmic Web. In the present analysis, we 
computed three hierarchical levels for practical reasons since 
further levels of substructure tend to be noisy due to reso- 
lution effects. We generate catalogues of walls and filaments 
and for all pixels in walls we also compute the distance to 
its closest adjacent filament. 

Figure [T] illustrates the close relation between LSS ge- 
ometry and local dynamics. The left panel shows the net- 
work of structures found with the SpineWeb along a thick 
slice. A small region containing a face-on wall at its centre is 
highlighted by a white square and the corresponding density 
field is shown in the central panel. Note how the substructure 
inside the wall contrasts with its adjacent empty voids. Even 
when the density of the wall is very low (8 ~ 0.5 — 1.0) it is 
denser than the interior of its adjacent voids which appear 
completely dark in the central panel. The right panel shows 
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Figure 1. Local dynamics inside a cosmological wall. The three panels show different aspects of the geometry and dynamics of the 
Cosmic Web. Left panel: network of walls, filaments and voids from the SpineWeb method shown in a thick slice of 1.36 Mpc/fc. Central 
panel: Density field inside the zoomed region highlighted in the left figure. Right panel: local dynamics from the residual velocity field at 
scales of 8 Mpc//i visualized with the LIC technique (see text for details). The lines represent streams following the velocity field. The 
black filled areas correspond to filamentary regions with overdensity 8 > 2. 



the dynamics inside the central wall. We visualize the veloc- 
ity field using the Line In tegral Convolution (LIC) technique 
(|Cabral and Leith II 19931 ) . The velocity field shown here cor- 
responds to the residual velocity field v" at a scale a defined 
as: vf — Vi—Vi* G(a), where Vi is the original velocity field 
decomposed in its three components and Vi * G(a) is its con- 
volution with a Gaussian kernel of scale a = 8 Mpc//i. The 
velocity flows delineated by the LIC streams closely resem- 
ble the inverse dynamical description given in section 12.11 
matter flows out of walls into their adjacent filaments and 
from there to the clusters. The outflow pattern at the inte- 
rior of walls is remarkably coherent and offers a geometric 
origin for an extra expansion factor. This outflow is sim- 
ply the result of the LSS geometry-dynamics as such it is 
independent of the background cosmology. 



6 THE HUBBLE FLOW IN WALLS 

Having defined the halo catalog we proceed to select those 
halos inside walls and filaments. Central MW halos in walls 
were restricted to be at least 3 Mpc away from their clos- 
est filament. No restriction was imposed on the neighbour 
haloes, except i> max > 20 km s _1 /i _1 . We also computed 
the local density fluctuation (S = pjp) around each halo 
from the dark matter density distribution using a Gaussian 
kernel with oq — 3 Mpc. We defined low (S < 1.3) and 
high (S > 2) density environments roughly corresponding 
to voids-walls and filaments-clusters respectively (see Fig. 
[1]). For each MW halo we compute line-of-sight velocities as 
v loB — r . iiMw-noi, where r is the unit vector between the 
central MW halo and its neighbour halo and e M w-n»i is the 
relative velocity between the two haloes. The radial veloc- 
ity is then used to compute the distance-velocity diagrams 
described in the following sections. 



6.1 Cold and super-cold Hubble flows in walls 

The Hubble flow for haloes in walls is shown in the distance- 
velocity diagrams of Fig. [2] for three increasingly cold and 



super-Hubble cases (compare to Fig. [3]). The flow for haloes 
in walls (left panel) has a mean velocity dispersion of an — 
55 km s _1 in the range 1 < r < 3 Mpc after which point 
we observe a broadening in the dispersion up to an = 180 
km s _1 at r = 5 Mpc. The ratio between the simulated and 
observed mean dispersions is a Bilri /a abls — 1.63. When the 
low-density condition is included (central panel) the Hubble 
flow becomes significantly colder with o~h — 30 km s _1 at 
1 < r < 3 Mpc. The broadening in an at larger distances 
is also smaller with an = 114 km s _1 at r = 5 Mpc and 
the ratio between mean dispersions is a" sim /a- ob5 = 1.07. The 
agreement of this particular sample and the observed galaxy 
population is remarkable with very similar dispersion across 
the full range of the diagram and even showing the same 
broadening at r ~ 3 Mpc. Finally, if one further restricts the 
neighbour haloes to be located in the same parent wall as the 
central MW halo (right panel) we then have the extreme case 
of super-cold flows with mean velocity dispersion oh — 25 
km s" 1 in the range 1 < r < 3 Mpc. The low velocity 
dispersion is observed even at the boundaries of the largest 
walls in our sample with an = 50.9 km s _1 at r = 5 Mpc. 
The dispersion in this sample is significantly lower than the 
observed by almost a factor of two (a Birn /a obB = 0.63). 



6.2 Super Hubble flows in walls 

The internal expansion of the walls is also observed in the 
distance velocity diagrams of Fig. [2] The flow for haloes in 
all cases is sub-Hubble at small scales r < 2 Mpc reflect- 
ing the decelerating influence of the central halo. For haloes 
in walls (left panel) the flow is marginally super-Hubble at 
r ~ 2.5 Mpc. Once we include the low-density condition 
(central panel) the super-Hubble character of the flow be- 
comes more clear. At distances r > 2 Mpc we have a local 
Hubble constant of h toca i = 0.77. When we consider only 
neighbour haloes inside the same parent wall (right panel), 
then the effect of the local LSS geometry-dynamics domi- 
nates the distance-velocity diagram. The local Hubble con- 
stant is hioad = 1-13 at scales r > 2 Mpc, which is 150% 
larger than the global Hubble constant given by the back- 
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Cold flow + marginal super Hubble 



Cold flow + super Hubble 
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Figure 2. Cold and super-Huble flows inside walls. We show the distance-velocity diagram for central "Milky Way" haloes inside walls 
considering three increasingly "cold" neighbourhoods: Top left panel shows the contribution from all haloes around central halos in walls. 
Top central panel is the same as left panel but considering only central halos with 8 < 1.3. Top right panel shows the case of 5 < 1.3 
and including only the contribution from neighbour haloes inside the same parent wall as the centre halo. The distribution of neighbour 
haloes is binned and presented as grayscale in the background. Blue line corresponds to the mean in the corresponding bin computed 
from a Gaussian fit. Blue dashed lines ind icate the standard deviation computed from the Gaussian fit. Filled circles correspond to 
galaxies given in lKarachentsev et al.l l|2009l 'l for r < 3 Mpc and Karachcn tsev et al.l l|2003h for r > 3 Mpc. The white solid line indicates 
a pure Hubble flow with h = 0.73. The radius of zero velocity Ro is a l so shown as vertical black dotted lines. Bottom panels show the 
corresponding ratios between the dispersion from N-body (cr a i m ) and the observed one (cr oba ). Values larger than one (hot) are plotted in 
red and blue (cold) otherwise. The blue color in the top panels indicates that for these curves <7 a i m /CT oba < 2 (see main text for details). 



ground cosmology (h — 0.73). This local Hubble constant is 
also significantly larger than the expected for g alaxies inside 
voids where h vaid = 1.2ft = 0.88 (|Schaadl20o"7f l. 



a factor of M oat /M rGa i = 337% compared to the low and 
high density samples (left and central panels of Fig. \3§ and 
M Bat /M re . a i = 462% compared to haloes in walls. 



6.3 Radius of zero velocity (Ro) 

We explored the effect of the local geometry-dynamics on 
the radius at which no recession velocity is measured (Ro)- 
This radius depends to first order on the mass of the central 
halo as Ro oc M 1 ^ 3 . The value of Ro is basically indepen- 
dent of the local density (see Fig. [3] left and central panels) 
at Ro — 1 Mpc. However, when one divides the MW halo 
sample in terms of their surrounding LSS geometry, the ef- 
fect of the local LSS dynamics emerges. For MW halos in 
walls Ro ~ 0.9 Mpc (see Fig. 0. This reflects the fact that 
neighbour haloes recede faster than the global Hubble flow, 
thus making the central halo appear less massive. The differ- 
ence in mass from no LSS-geometry distinction and haloes 
in walls is M wall /M all oc J?o,wau/-Ro >aJ i - l 3 /0-9 3 = 0.73. 
Without considering the local geometry for central halos in 
walls, one can underestimate its mass by as much as ~ 27%. 

An extreme example of the effect of local geometry- 
dynamics on the value of Ro is shown in the right panel of 
Fig. [3] where we computed the distance- velocity diagram for 
central haloes located in filaments and high density regions 
and considering only neighbouring haloes located in walls. 
This case is the opposite dynamically to the one of haloes 
in walls. The infall of haloes in filaments along the adjacent 
walls is reflected in the extremely high value of Ro- This 
effect is independent of the mass of the central halo (compare 
to the central panel of Fig. [3]) and if used to estimate the 
mass of the central MW halo can overestimate its mass by 



6.4 Effect of density on the Hubble flow 

Figure [3] shows the effect of local density on the distance- 
velocity diagram. Low density regions (left panel) have low 
velocity dispersions with <th = 63.2 km s^ 1 in the range 
1 < r < 3 Mpc. The value of an remains relatively constant 
for all distances. The mean velocity converges to the global 
Hubble flow at r ~ 4 Mpc. High density regions are sub- 
Hubble and significantly hotter with an = 153.8 km s _1 
going as high as an = 212.1 km/s at r = 5 Mpc. 



7 CONCLUSIONS AND FUTURE WORK 

The coldness of the local Hubble flow is a longstanding prob- 
lem in ACDM cosmology. We have shown that by taking ac- 
count of the void/wall/filament nature of large-scale struc- 
ture, one can naturally reproduce the observed low disper- 
sion in the local Hubble flow for galaxies located in local- 
group LSS configurations, corresponding to walls. The dis- 
persion of the local Hubble flow (1 < 7? < 3 Mpc) around 
simulated "Milky Way" haloes which are located at the cen- 
tre of low-density cosmological walls is found to be o~h ~30 
km s _1 (Fig. [2l central panel), in excellent agreement with 
observations. The expansion inside walls is also reflected in 
the enhanced value of the measured local Hubble constant. 
For "Milky Way" haloes inside walls, we find super-Hubble 



Cold flow 



Locally Cold Flows from Large-Scale Structure 5 

Hot flow + sub-Hubble Hot flow + sub-Hubble + High Ro 




Figure 3. Distance-velocity diagram for MW haloes for three different environments: left panels: central haloes in low density regions. 
Central panels: central haloes in high density regions. Right panels: central haloes in high density regions and located inside filaments 
and considering only neighbour haloes located inside walls. The figure description is the same as in Fig. [3] The blue color in the top-left 
panel indicates that for this curve cr siln /(j obs < 2. Top central and top right panels are colored in red indicating <T sim /(T b s > 2. 



flows with ftiocai — 0.77 — 1.13. Ignoring the dynamical ef- 
fects of the local wall, the mass of the Milky Way estimated 
from Ro can be underestimated by as much as ~ 30%. In a 
future paper we will extend our study to include voids and 
filaments. 

The strong constraint on the LSS dynamics imposed 
by the complementary LSS geometry can also be exploited 
to improve current techniques of velocity field reconstruc- 
tion. A full LSS morphological decomposition of the local 
supercluster can allow us to obtain a more complete char- 
acterization of our local velocity field. 
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